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ABSTRACT 

The 3D Radioastronomical Approach to Doppler tomography has been used to study the Ha emis- 
sion sources in U Coronae Borealis. These 3D tomograms provide greater resolution than the projected 
2D version and highhght the jet-Hke gas flows in the 14 direction transverse to the orbital plane. In 
this paper, the 3D tomograms are compared at two distinct epochs when U CrB was in the disk-like 
state (1993 data) and the strcam-likc state (1994 data). Both states display a prominent emission 
source, the circumprimary bulge, which is produced when the gas stream strikes the photosphere of 
the mass-gainer. This source is detected within = ±150 km s~^, and demonstrates that the bulge 
is not confined to the orbital plane although it achieves maximum strength near 1^=0 km s~^. Other 
emission sources include the stream-star and stream-disk shocks and a Localized Region (LR) where 
the circling disk material strikes the incoming gas stream. The LR has Vz velocities of 200 to 500 km 
s~^ in the disk-like state. The disk emission is seen over a range of Vz velocities, and there is evidence 
that the disk is inclined to the orbital plane or may have two arms. The gas stream flows along its 
predicted trajectory in the stream-hke state, and a comparison with the disk-like state suggests that 
the gas stream has a higher density than the disk in both states of this binary. 

Subject headings: techniques: image processing accretion, accretion disks - stars: binaries: close - 
(stars:) binaries: eclipsing (stars:) circumstellar matter - stars: imaging - stars: 
individual (U Coronae Borealis) 



1. INTRODUCTION 

The two-dimensional Doppler tomograms of the Algol- 
type binaries (Richards. Albright & Bowles 1995; Al- 
bright & Richards 1996; Richards 2004) display a di- 
verse range of circumstellar structures. These include a 
gas stream, accretion annulus, transient accretion disk, 
shock regions and sometimes a chromospheric emission 
source in the short-period (P < 5 days) Algols, as well 
as a classical accretion disk in the long-period (P > 5 - 6 
days) Algols. Some systems, like U Coronae Borealis (U 
CrB) and U Sagittae (U Sge), alternate between stream- 
like and disk-like states and display changes that can 
occur within a year, and sometimes overnight (Richards 
& Albright 1999). These 2D Doppler tomograms are 
created by assuming that only one V, velocity value is 
represented, namely, Vz = km s~^. However, since the 
2D map is a superposition of all the slices, then all of the 
Vz velocity information is embedded in this integral 2D 
map. The primary aim of this work is to demonstrate 
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that 3D Doppler tomograms can be generated by assum- 
ing a grid of Vz values and that the resulting 3D maps 
are consistent with the 2D Doppler tomograms. 

The short-period Algols represent the class of direct 

impact binaries in which the gas stream from the Li point 
strikes the photosphere of the mass gaining star directly. 
This occurs because the radius of the mass gainer is large 
relative to the binary separation, so the gas stream is 
not deflected sufficiently by the Coriolis force to avoid 
c(mtact with the star (Richards & Albright 1999). Two- 
dimensional hydrodynamic simulations of direct impact 
binaries, like U CrB, have illustrated how the high speed 
(~ 500 — 600 km s~^) gas stream is shocked and heated 
by the impact with the slowly rotating (?;sini < 150 
km s~^) mass gainer (Blondin, Richards & Malinkowski 
1995; Richards & Rathff 1998). Some of the gas from 
the impact/shock region is released into orbits that may 
repeatedly splash off the star until energy losses due to 
these collisions cause the gas to be accreted onto the 
star (Smak 1989). The dynamics of the impact (Ul- 
rich & Burger 1976; Richards 1992) suggest that the 
shock should generate gas motions transverse to the or- 
bital plane. Richards & Ratliff (1998) simulated both 
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the Cartesian and velocity evolution of the gas distri- 
bution in direct impact systems and created synthetic 
2D Doppler tomograms to compare with the 2D tomo- 
grams derived from observed spectra. They found that 
the change between the disk-like and stream-like states 
could be explained by a change in the mass transfer rate 
in the binary. Three-dimensional hydrodynamic simu- 
lations of the (3 Lyr eclipsing binary by Bisikalo et al. 
(2000b) found jets moving perpendicular to the orbital 
plane. So, substantial disk structure transverse to the 
orbital plane can form naturally in these binaries. 

The extension to 3D Doppler tomography was first for- 
mulated in the classical paper by Marsh & Horne (1988), 
and more recently discussed briefly by Steeghs (2004) 
and Marsh (2005), who suggested that "out-of-plane 
motions" needed to be included to properly explain the 
observed spectra. However, the extension to the Vz di- 
rection remained unexplored until the work of Agafonov, 
Richards & Sharova (2006), which outlined a new pro- 
cedure for creating 3D Doppler tomograms. The first 
3D Doppler tomogram of a binary was constructed by 
Agafonov, Richards & Sharova (2006) from forty-seven 
Ha spectra of the Algol-type binary U CrB. They used 
the reconstruction method called the Radioastronomi- 
cal Approach (RA), which was developed by Agafonov 
(2004a, b); Agafonov & Sharova (2005a, b) for the case 
when there is a small number of projections (called few 
projections tomography; see Section 3.1). The 3D to- 
mography of U CrB displayed an enhanced image of the 
gas stream along the ballistic trajectory and an equato- 
rial emission source centered on the primary star that 
were more distinct than the image produced by standard 
2D Doppler tomography (Richards, Albright & Bowles 
1995; Richards 2001). A high-velocity gas jet with a 
high emission intensity was also discovered moving in 
the direction transverse to the orbital plane of the sys- 
tem with a high Vz velocity component. This jet could 
not have been detected previously in the standard 2D 
tomograms because all velocity information is projected 
onto the (T4,l^j/) plane only, and hence the 14 compo- 
nents cannot be represented in a 2D map. 

In Cartesian coordinates, a "jet" is an elongated fea- 
ture formed by gas motions within the structure. In this 
work, a jet is a luminous (high intensity) feature with 
the following properties: (1) It can be elongated in the 
(14;,Vy) velocity planes that have Vz components close to 
zero. (2) The T4 and/or Vy velocity values of the fea- 
ture should be significantly different from zero, since all 
(14,Vy) planes (for all values of Vz) are associated with 
the rotation of the binary. (3) At the same time, the jet 
feature might not be elongated in the {Vx,Vy,Vz) veloc- 
ity space if it has a high Vz velocity component since this 
feature would move in a direction transverse (or with sig- 
nificant inclination) to the direction of the orbital plane 
of the system. So, it is reasonable to call the gas stream 
from the LI point a jet even if it has small Vz velocities. 
Other jet-like features could have velocities in the Vz di- 
rection. In general, any gas flows in the {Vx,Vy) plane at 
Vz=0 should move in planes that are collinear with the 
{x,y) orbital plane. Moreover, since gravitational fields 
typically dominate the gas flows, the gas velocities in the 
iVx,Vy) plane at Vz=0 should naturally coincide with the 
Cartesian (x,y) orbital plane of the binary. So, reason- 



able arguments can be made about the nature of the flow 
in the orbital plane based on the flow behavior in the 
"central velocity plane" deflned by Vz=0, even though 
the velocity fields of the composite emission sources are 

not well-known. 

In this paper, we have continued our effort to under- 
stand the 3D structure of accretion fiows by: 1) creat- 
ing a procedure for constructing the 3D Doppler tomo- 
gram from the ID projections, 2) constructing the 3D 
tomograms and identifying the optimum way to display 
these tomograms, 3) identifying the features in the 3D 
images and relating them to those seen in the 2D tomo- 
grams, and 4) providing an initial physical interpretation 
of these features. In addition, we have used 3D tomog- 
raphy to study the variability of the gas flows of U CrB 
in its disk-like and stream-like states. A lot of attention 
was focussed on the fine structure extraction of the to- 
mogram and the analysis of the features of U CrB when 
it was in the disk-like state. The visualization and anal- 
ysis used in the 3D Doppler imaging process have been 
advanced in this work by using the RA reconstruction 
procedure. It is now possible with our improved under- 
standing of the 3D velocity images to make preliminary 
interpretations of the 3D accretion flows. The study of 
3D Doppler tomograms of several binaries in a variety of 
states will help us to gain a better understanding of the 
significance of these transverse flows and their role in the 
accretion process. 

2. 2D DOPPLER TOMOGRAPHY OF U CRB 

The Algol-type binary U Coronae Borealis is an inter- 
acting system in which the cooler GO III-IV secondary 
star is transferring mass to its B6 main-sequence com- 
panion (the primary) through Roche lobe overflow. U 
CrB is a short-period Algol (P=3.4522 days) with an or- 
bital inclination of 79.1° (Cester et al. 1977) and mass 
ratio g=0.29 (Batten & Tomkin 1981). The masses of 
the primary and the secondary are equal to 4.8 and 1.40 
Mq and their radii are equal to 3.0 and 4.6 Rq, respec- 
tively. The systemic velocity of the system is Vo= -6.7 
km and the velocity semi- amplitude of the primary 
is Kp = 59.7 km s~^. One hundred and thirty-nine Ha 
(6562.8 A) spectra with a resolution of 0.093 A/pixel 
or 4.3 km s~^ were collected in 1993 with the NSO - 
McMath-Pierce telescope, and forty-seven spectra with a 
resolution of 0.166 A/pixel or 7.6 km s^^ were collected 
in 1994 with the KPNO 0.9m Coude Feed Telescope (see 
Richards & Albright 1999). The spectra were collected 
at closely spaced positions around the entire orbit of the 
binary. Since the observed spectra are dominated by the 
spectrum of the B6 primary star, the spectral contri- 
butions of the stars were subtracted from the observed 
spectra to create difference spectra (see Richards 2004 
for details). In this procedure, model atmospheres calcu- 
lations were used to represent the photospheric contribu- 
tions of the stars. The difference spectra, which display 
enhanced emission and absorption contributions from the 
circumstellar gas, were processed for image reconstruc- 
tion. 

The earlier standard 2D Doppler tomograms of U CrB 
1993 (in the disk-Uke state) and U CrB 1994 (in the 
stream-like state) taken from Richards (2001) are shown 
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in Figure 1. These 2D images were compared with the 
3D Doppler tomograms based on the same data. The 
1993 tomogram shows the presence of the following fea- 
tures (Richards 2001): an accretion disk around the 
primary star, the gas stream, the stream-star impact re- 
gion, a circumprimary emission source, and an absorp- 
tion zone (region dominated by absorption). However, 
Figure 1 shows that the contrast of these separate fea- 
tures is not high because these sources are superimposed 
on each other in the 2D velocity image. The standard 
2D Doppler tomogram reflects the summed emission on 
the {Vx,Vy) plane and is somewhat distorted by dispersed 
rays of weak emission, which are a direct consequence of 
the filtered back-projection method. The second frame of 
Figure 1 shows the gas stream along the predicted grav- 
itational trajectory as well as circumprimary emission. 
However, these features are merged instead of appearing 
as distinct sources. 

3. 3D DOPPLER IMAGING 

The reconstruction of the 3D Doppler tomogram of U 
Corouae Borealis based on the 1993 Ha line profiles as 
well as the identification of emission features are pre- 
sented in this section. This information has been used to 
provide new information about the binary. 

3.1. Reconstruction and Visualization of the 3D 
Tomogram 

The Filtered Back-Projection Method applied in com- 
puted tomography was first described in the classical pa- 
per by Bracewell & Riddle (1967). It had wide use in 
the solution of a variety of tomography problems and 
it was also applied to astronomical applications. It is 
now widely used in Doppler tomography for image re- 
construction. Many Doppler maps have been created us- 
ing the Maximum Entropy Method (MEM) in which the 
deconvolved image is passed through a series of statisti- 
cal tests to ensure that it is consistent with the mea- 
sured data (e.g., Marsh & Home 1988). The MEM 
then selects one image from the most feasible set of im- 
ages that describe the data equally well. However, the 
simpler back-projection method has been used to create 
Doppler images of cataclysmic variables and Algol bi- 
naries (e.g., Kaitchuck et al. 1994; Richards, Albright 
& Bowles 1995, and others). The standard version of 
this back-projection technique maps the Vx and Vy com- 
ponents only and produces images in 2D velocity space. 
However, the gas fiows in the directions away from the 
orbital plane may be significant for many binaries. The 
development of the 3D version of Doppler tomography 
has now permitted the study of accretion structures be- 
yond the orbital plane by incorporating the Vz velocity 
component of the moving gas flows. 

The 3D reconstruction applied in this work is called the 
Radioastronomical Approach (RA). It was developed by 
Agafonov & Sharova (2005a,b). The theoretical possi- 
bility of image reconstruction from projections was de- 
scribed by J. Radon in 1917 (Radon 1917). However, 
the mathematical equations of the Radon transform and 
its inverse are formulated for the ideal case with infl- 
nite numbers of projections and inflnitely high resolu- 
tion, covering the entire range of the frequency domain. 



The widely-used Filtered Back Projection (FBP) method 
is based on the work of Bracewell & Riddle (1967) 
and permits the application of a flnite number of pro- 
jections; with the image reconstruction occurring over 
a wide range of spatial frequencies in the 2D Fourier 
domain. Equations in Bracewell & Riddle (1967) can 
be used to estimate the number of projections required 
to produce the best reconstructed image and includes 
the filtering of redundant high spatial frequencies. The 
more general Radioastronomical Approach is based on 
a different set of initial conditions than FBP. It per- 
mits 3D tomographic reconstruction (see Agafonov & 
Sharova 2005b for details). This method is especially 
useful when the number of projections is relatively small 
because it can achieve the same resolution as FBP while 
using fewer profiles. The RA can be appHed in the usual 
cartesian coordinate space to solve a variety of problems. 
So, in the case of Doppler tomography the main distin- 
guishing features of the RA procedure are: (1) the so- 
lution of the inverse problem in 3D velocity space, (2) 
the introduction of the summarized image and summa- 
rized point spread function (SPSF), and (3) the decon- 
volution of the image using the CLEAN technique to re- 
move the responses from the sidelobes of the SPSF. The 
summarized 3D image and summarized 3D point spread 
function were calculated based on the line profiles and 
one-dimensional (ID) transfer functions by incorporat- 
ing the angles corresponding to the orbital inclination 
and orbital phases. The required version of the 3D recon- 
struction corresponds to the 3D id variant described by 
Agafonov & Sharova (200-5a) and Agafonov & Sharova 
(2005b). The SPSF was calculated on the set of knife- 
beam ID transfer functions, which determine the velocity 
resolution of the profiles. This 3Di]3 notation represents 
the 3D reconstruction based on the ID profiles. Here, the 
knife-beam function is the function or beam that trans- 
fers the intensity from the 3D figure to the ID profile or 
backwards from the ID to the 3D image (see Agafonov 
& Sharova 2005a). Thus, the RA reconstruction is sig- 
nificantly different from all other reconstruction schemes 
associated with computed tomography and currently ap- 
plied to medicine, since those techniques derive the 3D 
image from the set of 2D slices and the reconstruction of 
these 2D slices is based on the FBP method (e.g., Kalen- 
der 2005). 

In the reconstruction of the 3D Doppler tomograms, 

we analyzed the same set of Ha line profiles that were 
used by Richards (2001) to construct the standard 2D 
Doppler tomogram. The 3D Doppler tomogram was cal- 
culated as a set of numerical values in the cells of a cube, 
with dimensions {Vx,Vy,Vz) having values ranging from 
-700 to +700 km s^^. The maximum of this function 
was normalized to the unit to permit the visualization of 
the image based on any slice taken from any direction. 
To simplify the physical situation, the {Vx,Vy) slices are 
displayed in the horizontal plane (to be consistent with 
the direction of the orbital plane of the binary) , and the 
slices containing the axis are in the vertical direction. 

Figure 2 displays 15 slices in the {Vx,Vy) plane for val- 
ues of Vz ranging from —420 km s~^ to +420 km s~^, in 
60 km s~^ steps. The superposition of the color and con- 
tour images allows us to emphasize the main features in 
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each Vz sUce of the 3D image. Figure 3 displays enlarged 
images of the four most interesting Vz slices. Other slices 
at positive Vz values are shown in Figures 4 and 5, and 

several important cross-sections in the (VyjVz) plane are 
presented in Figures 6 and 7. Figures 5 and 7 are in- 
cluded here to show the 3D slices produced from the 
1994 spectra by Agafonov, Richards & Sharova (2006). 
Figure 8 illustrates the comparison between the distinct 
disk-like and stream-like distributions derived from the 
1993 and 1994 spectra. In Figures 2-8, the blue contours 
represent absorption while red contours represent emis- 
sion. Finally, Figure 9 represents the proposed model of 
the U CrB binary. 

3.2. Description and Interpretation of Emission 
Features 

Both Figures 2 and 3 demonstrate that the most promi- 
nent emission in the UCrB (1993) image arises from the 
accretion disk around the primary star, as expected from 
the 2D image. The disk emission strengthens as we move 
from Vz= —420 km towards the Vz— km s^^ slice 
and then weakens again hy Vz= +420 km s~^. This Vz 
velocity range will decrease with the adjustment for the 
lower resolution in the Vz direction. The disk emission 
has an asymmetrical structure, as though composed of 
two parts with different velocities. Figure 2 shows that 
the maximum intensity of the mass gaining star occurs 
near the central Vz velocities (near 1^=0). However, the 
maximum emission intensity in the 3D tomogram {Imax 
= 1.0) was found in the Vz= —20 km s^^ slice (Figure 3), 
and corresponds to the velocity location of the primary 
star. The intensities of the mass loser and mass gainer 
are aligned in a symmetrical way along the 14=0 line. 
The intensities of these sources was strongest near Vz=0, 
decreased in both directions towards Vz= ±180 km s^^, 
and was negligible at higher absolute Vz velocities. 

We consider ten bright emission features that are dis- 
tinctly visible in the tomograms. Most of these features 
were identified earlier in the 2D Doppler tomograms, 
and from comparisons between the predicted locations 
of these structures in both velocity-space and Cartesian 
space. The ten features are outlined in Table 1 along with 
their typical velocity ranges. Eight of these features are 
identified in Figure 3, which shows the four most inter- 
esting slices corresponding to Vz = —240,-20,240 and 
420 km s~^ in the 1993 tomogram. The features are (1) 
circumprimary emission and accretion annulus, (2) chro- 
mospheric emission associated with the of cool secondary 
donor star, (3) the stream-disk impact region (hot-spot) 
in the middle of the predicted path of the stream, (4) the 
stream-star impact region, where the gas stream strikes 
the stellar surface, (5) the locus of a Keplerian accre- 
tion disk, (6) a bright feature near the center of the to- 
mogram, previously identified as the L(K:alized region by 
Richards (1992), (7) another feature which may be asso- 
ciated with the Localized Region, and (8) a feature near 
the velocity of the mass loser, but with high Vz velocities. 
The sources labeled 3, 6, 7, 8 dominate in brightness at 
the higher positive Vz velocities. These sources are illus- 
trated in Figure 4, for slices from Vz =240 to 540 km s~^. 
Two additional sources were identified in the 1994 tomo- 
gram of Agafonov, Richards & Sharova (2006) and are 



shown in Figure 5: (9) the gas stream fiowing along the 
predicted ballistic trajectory, and (10) a high velocity 
stream (or jet) moving in the transverse direction rela- 
tive to the orbital plane of the binary. 

The velocity resolutions in the {Vz,Vx) plane and the 
{Vz,Vy) plane depend on the orbital inclination of the bi- 
nary and are equal only for an inclination of 45°. For 
high values of inclination, the resolution in the Vz di- 
rection is not as good as in the other dimensions, and 
gets worse as the values of Vz increase. Since the or- 
bital inclination of U CrB is high (i = 79.1°), the 3D 
tomogram was restored with a resolution of 20 km s^^ in 
the Vx and Vy directions and 70 km s~^ in Vz direction. 
So, the features on the tomogram look elongated in the 
Vz direction because the resolution in the Vz direction 
is approximately 3.5 times lower than in the V^ and Vy 
directions (see Figures 6 and 7). The full dimensions of 
all features are given in Table 1 and throughout the text. 

The possibility of an adjustment for the lower resolu- 
tion in the Vz direction depends on whether a source is 
compact or elongated. A source is considered to be com- 
pact if it has one expressed intensity maximum and a 
narrow range of velocities, while an extended source is 
characterized by a wide range of velocities in more than 
one direction {Vx,Vy or Vz)- For example, the circumpri- 
mary emission (Feature 1) is obviously compact because 
it has only one intensity maximum and a narrow range 
of velocities. Its isolines (equal intensity lines) have cir- 
cular shapes because of equal velocities in the V^ and Vy 
directions (Figure 2 or 3). Feature 1 is artificially elon- 
gated in the Vz direction (in Figure 7, right frame) only 
because of the poorer resolution in that dimension. So, 
the measured Vz velocities for compact sources should 
be scaled by the factor of 1/3.5 (= 0.29) before direct 
comparison with the Vx or Vy velocities. The deconvolu- 
tion procedure is more complicated for extended sources 
(e.g.. Features 6 - 8), so any features with the ratio Vz/Vx 
or Vz/Vy greater than 3.5 are considered to elongated in 
the Vz direction. Consequently, the shapes of the isoUnes 
in the {Vy,Vz) or {Vx,Vz) planes can be compared with 
those in the {Vx,Vy) plane once we take into account the 
factor of 0.29 scaling of the artificially stretched shapes 
of all the components in the Vz direction. 

The compact blue absorption feature seen in the lower 
left quadrant of Figure 2 from Vz = -180 to -1-420 km s~^ 

(see also Figure 8, lower right frame) was called the ab- 
sorption zone by Richards (2001, 2004), and is associated 
with the locus where the gas temperature becomes too 
high to emit at Ha. A structure at a similar location 
was identified in the ultraviolet spectrum of U Sge by 
Kempner & Richards (1999). 

A short description of each feature is provided below. 
Figures 2-7 have been used to examine and interpret 
the features. 

1.) The circumprimary emission (centered on the mass 
gainer) and an accretion annulus (a circular region of low 
velocity) are labeled (1) in Figure 3. These structures are 
centered along the Vy axis {Vx = km s^^) and slightly 
below the Vx axis in all the slices from Vz= —120 to -1-120 
km s~^. The full extent of this structure is approximately 
360 km s~^ in the Vz direction, and ranges from about 
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T4= -180 to +180 km s"^ (see Figure 2), while its width 
is only 120 km s~^ in the {Vx,Vy) plane. Since the reso- 
lution in the T4 direction is 3.5 times worse than that in 
the {Vx,Vy) plane (Section 3.1), the corrected Vz width 
should be (360 km s-i/3.5) -100 km s'^ (or 14=-50 to 
Vz= +50 km s~^), which is comparable to the width of 
the structure in the {Vx,Vy) plane. This result suggests 
that the circumprimaxy structure is nearly spherically 
symmetric in the velocity domain. 

2. ) Chromo spheric emission from the cool donor star 
(the secondary) is shown at (2). Its Vz velocity ranges 
from Vz= —120 to +120 km s~ , as shown in Figure 2. 

3. ) A bright (hot) spot in the middle of the gas stream 
is located at (3) and represents the site where the gas 
stream makes impact with the outer edge of the accre- 
tion disk. It is detected in the slices from Vz— +180 to 
+540 km s~^. However, emission along the predicted 
gravitational trajectory of the gas stream looks faint in 
all slices based on 1993 spectra. This stream emission 
is clearly brighter in the 1994 tomogram (see Figure 5). 
Feature 3 may also represent gas ejected out of the or- 
bital plane by the impact with the disk. This ejected gas 
has high positive Vz velocities from 240 to 540 km s^^, 
centered on 420 km s^^, near Vx = -360 km s^^ (see Fig- 
ure 6 and the slice for Vz=420 km s^^ in Figure 3). It 
is apparent that the ejected gas then flowed within the 
predicted locus of the accretion disk, along the trajectory 
marked in Figure 3 (bold dashed line), and with a lower 
Vz velocity of 240 km s~^. As this gas blends with the 
disk (Feature 5), it acquires negative Vz velocities up to 
-200 km s~^. This transition is evident in Figure 3 (slice 
Vz= -240 km s-i). 

4. ) The location of the stream-star impact region, where 
the gas stream strikes the stellar surface, is found at (4) 
near the location where the gas stream trajectory inter- 
sects the inner part of the disk. Its velocity ranges from 

T4= -300 to +60 km s'^. 

5. ) Beyond the impact/splash zone, the gas motion at 
(5) is consistent with the predicted locus of the accretion 
disk and rotates like a disk. Its Vz velocities ranged from 
-420 to +420 km s~^ (see Figures 2 - 4), so this disk has 
substantial flow in the vertical direction, although the ve- 
locity range will be slightly smaller after the adjustment 
for the lower resolution in the Vz direction. This result 
was expected from studies of spectra collected during pri- 
mary eclipse. This disk may be turbulent with vortical 
movements, especially if it is inclined. We can identify 
the stream-disk and stream star regions even when the 
disk emission is prominent if the disk has a lower den- 
sity than these regions or if it is inclined to the orbital 
plane of the binary. Some evidence in favor of an in- 
clined disk is provided by Figure 2, in which the lower 
left part of the disk is visible mainly for negative Vz slices 
and the emission seems to shift to the lower right part 
of the disk for positive Vz slices. In addition, one part 
of the stream (flowing along the gravitational trajectory) 
makes impact with the disk, and subsequently, the bright 
feature 3 (hot spot along the stream) has positive veloc- 
ities (so it is moving away from us). At the same time, 
the stream-star impact region (Feature 4) has negative 
Vz velocities, so this transverse region is moving towards 
us (Figure 6) . The emission intensity of the disk is not as 



strong as some of the other features, but it is evident that 
the gas flow in the disk spreads around the primary star. 
Moreover, we cannot exclude the existence of two arms 

within the disk with different sets of velocities, as was 
predicted by Sawada et al. (1986) and Spruit (1987) 
and then observed in the case of the cataclysmic variable 
IP Peg by Steeghs et al. (1997) and others. 

6 & 7.) The two concentrated regions labeled (6) and (7) 
are found near the center of the tomogram in the same 

general region that was identified as the Localized Re- 
gion in Cartesian models of (3 Per (Richards 1992, 1993), 
RW Tau (Vesper & Honeycutt 1993), and other direct- 
impact systems. These regions are detected at high pos- 
itive Vz velocities ranging from +240 to +540 km s~^. 
At the highest Vz velocity, they become part of the same 
structure. As explained by Richards (1992), this overall 
structure is created by the natural flow within the disk 
when it circles the primary and strikes the incoming gas 
stream, to produce a disk-stream shock region. A simi- 
lar structure consisting of jets moving perpendicularly to 
the orbital plane was also identifled from hydrodynamic 
models of the eclipsing binary, (3 Lyr A (Bisikalo et al. 
2000b). These regions can be examined in greater in de- 
tail in Figure 4, which shows fragments of six slices cor- 
responding to high positive Vz up to +540 km s~^. They 
can also be seen in Figure 6 which shows the (K/,Vz) 
cross-sections. This flgure shows that features 6 and 

7 have approximately the same Vz velocities {Vz ~400 
km s~^) as feature 3 (which represents the stream-disk 
impact in the middle of the stream). Feature 3 is seen 
in slices Vx= -360 and -300 km s~ while Feature 6 is 
shown in slices Vx= -20 and km s^^. Although Fea- 
tures 3 and 6 have similar Vz velocities, they have dif- 
ferent Vx velocities; the gas velocity vector for Feature 
6 is almost perpendicular to the central velocity plane 
(and hence the orbital plane), while Feature 3 has two 
velocity components {Vx and Vy) similar to the disk, its 
Vz velocity component has a high value of ~400 km s~^. 
As a result, it is reasonable to conclude that Feature 3 is 
the part of the disk gas that was ejected out of the or- 
bital plane after the impact with the stream. The regions 
6 and 7 could explain why the disk gas breaks up into 
clumps and changes its velocity direction near the site 
of the Localized region where the circling disk gas makes 
impact with the trajectory of the gas stream (Richards 
1992; Kuznetsov et al. 2001; Bisikalo et al. 2000a). 

8.) The feature labeled (8) has a velocity similar to the 

Vx and Vy velocities of the donor star, but with high pos- 
itive Vz velocities in the range from X4=+200 to +540 
km (sec Figures 2 6). This structure may be chro- 
mospheric in origin but the high velocities suggest that 
it could correspond to a jet or flare associated with the 
chromosphere of the secondary star. However, it is also 
likely that features 6, 7, and 8 may all be related to the 
interaction between the material in the disk and the in- 
ner edge of the stream. The Vz = 480 km s~^ slice in 
Figure 4 hints at a connection between regions 7 and 8, 
although they seem to be separate in other slices. The 
low velocity gas that circles the star may merge with the 
disk when it interacts with the gas stream (to become 
Features 6 & 7, and perhaps 8), or it could escape from 
the system in the Vz direction. 
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9. ) The feature labeled (9) represents the gas stream 
flowing along its predicted ballistic trajectory. It flows 
from the LI point towards the mass gaining star with 
speeds up to 450 km s^^. This feature is seen in Fig- 
ure 5 for the U CrB (1994) data (Agafonov, Richards & 
Sharova 2006). 

10. ) The feature labeled (10) is the high- velocity stream 
or jet. It is a compact source in the {Vx,Vy) plane but 
has a significant transverse flow in the {Vy,Vz) plane, with 
negative Vz velocities primarily from —400 to km s^^ 
(Agafonov, Richards & Sharova 2006). It can be seen in 
Figures 5 and 7 for the U CrB (1994) data. Gas speeds of 
up to 600 km s~^ were expected at the star stream impact 
site from dynamical arguments presented by Ulrich & 
Burger (1976); Richards (1992). 

3.3. Variability of the Gas Flows in U CrB 

The 3D Doppler tomogram of the U CrB binary based 
on the 1993 data (hereafter, 3D-1993) displays an ex- 
tensive accretion disk, a circumprimary emission source, 
a chromospheric emission source and a jet-hke structure 
associated with the secondary star, as well as stream- 
star and stream-disk shock regions. Albright & Richards 
(1996) suggested that the binary was in a disk-like state 
at this epoch. It is useful also to compare these results 
with those derived from the 1994 data (hereafter, 3D- 
1994). 

The 2D tomogram of U CrB based on 1994 data (Fig- 
ure 1) shows that the gas flow in the binary was in a 
stream-like state during that epoch (Richards, Albright 
& Bowles 1995). This 2D image is dominated by the gas 
stream which flows along the predicted gravitational tra- 
jectory from the LI point towards the mass-gaining star. 
Agafonov, Richards & Sharova (2006) used 3D Doppler 
tomography to produce a more distinct image of the gas 
stream than possible with the standard 2D method. The 
cascade visualization of the 3D- 1994 image is shown in 
Figure 5, which demonstrates that most of the gas stream 
flow (Feature 9) was in the central planes close toV^ = 0. 
It is obvious that the motion of any such gas must will be 
in the central velocity plane (and also the orbital plane) 
if it has no velocity component. A transverse flow of 
gas (a gas jet. Feature 10) was also found with a high 
Vz velocity component and two other weak components 
{Vx=-60 and Vy=0 km s^^; Figure 7, left). It is prob- 
able that a feature with such velocity components may 
have arisen from the expected location of the stream-star 
impact site (Agafonov, Richards & Sharova 2006). 

The cascade visualization of the 3D-1994 image is 
shown in Figure 5, which demonstrates that most of the 

gas stream flow (Feature 9) was in the central velocity 
plane defined by 14 = 0. It is obvious that the motion of 
any such gas must will be in the c;entral velocity plane if 
it has no Vz velocity component. A transverse fiow of gas 
(a gas jet. Feature 10) was also found with a high Vz ve- 
locity component and two other weak components (Ke= 
—60 km s~^ and Vy—0 km s"^; Figure 7, left). It is prob- 
able that a feature with such velocity components may 
have arisen from the expected location of the stream-star 
impact site (Agafonov, Richards & Sharova 2006). The 
{Vy,Vz) cross section of this gas jet is shown in Figure 



7 (left) for Vx= —60 km s~^, adapted from Agafonov, 
Richards & Sharova (2006). Figure 7 (right) also shows 
the {Vy,Vz) cross section (at T4= km s~^) of the emis- 
sion source which corresponds to an equatorial bulge on 
the primary star taken from the same 3D- 1994 image. 

The {Vy,Vz) cross sections for the 3D-1993 image (Fig- 
ure 6) can be compared with the 3D-1994 results shown 
in Figure 7. The center of the equatorial emission on the 
primary star (feature 1) is displaced at Vz = 80 km s~^ 
towards positive Vz velocities in 3D- 1994. but the same 
feature is located near the center in 3D- 1993. The dis- 
placement is only of -20 km s^^, which is not significant. 
However, the differences in the velocity structure of the 
transverse flows in the two states of the binary may be 
significant. 

The {Vy,Vz) cross section for 14 = —60 km s~^ for 
3D-1994 shows the transverse flow. The velocity of the 
transverse flow (Feature 10) in 1994 (Figure 7) ranges 
from Vz = —400 to km s^^, and the flow is nearly per- 
pendicular to the orbital plane. The gas stream impact 
site (Feature 4) in 3D- 1993 has a similar velocity range 
at the end of the stream in the slice for V^= -580 km s~^, 
but this gas does not flow in a perpendicular direction 
across the orbital plane. The correspondence between 
Features 4 and 10 is reasonable since this jet may have 
resulted from the splash at the impact site (see Agafonov, 
Richards & Sharova 2006). However, the other trans- 
verse gas flows for 3D- 1993 have the opposite direction 
relative to similar flows in 3D-1994. Feature 3 in the 
middle of the stream and features 6, 7, and 8 range from 
Vz = 200 to 500 km s~^ (see Figure 6). This is the main 
difference in the behavior of the transverse flows at the 
two epochs. 

The high velocity transverse jet (Feature 10; Figure 5) 
and the localized region (Feature 6; Figure 4) have simi- 
lar locations in the {Vx,Vy) plane, but the cross sectional 
views show that these features move in different direc- 
tions (see Figures 6 and 7) in the {Vy,Vz) plane. Feature 
6 has positive velocities with V^ 400 km s~^, and is 
flowing away from us, while Feature 10 has negative ve- 
locities with Vz ^ —200 km s^-', and is flowing towards 
us. However, it is interesting that Feature 10 (the high 
velocity transverse jet; Fig. 5b) and Feature 4 (the star- 
stream impact region; Fig. 5a) have the same velocity 
range from Vz = —400 to km s~^, which conflrms the 
conclusion by Agafonov, Richards & Sharova (2006) that 
the high velocity jet is located at the site of the impact 
between the gas stream and the star. 

It should also be noted that the range of Vz velocities 
for feature 9 (gas stream) is nearly symmetric about the 
Vz=0 direction in the 3D-1994 images (Figure 5). We 
also notice that the gas stream or jet (feature 10) flow- 
ing in Vz direction has negative Vz velocities, but later 
the velocities become positive suggesting motion of the 
gas around the binary. As a result, the circumprimary 
emission in 1994 has the displacement of the maximum 
intensity to T4=80 km s~^ (in positive values of Vz)- This 
fact was noted in the stream-state of the primary. 

While the image reconstruction process can now pro- 
duce 3D velocity images, we are still unable to create the 
desired 2D or 3D Cartesian images because the various 
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emission features have different velocity fields. However, 
an estimate of the size of the accretion disk can be made 
from information provided in Figure 1. Based on the 

point where the gas stream is truncated, we can pro- 
duce a Cartesian model to show the dominant emission 
sources in U CrB. The distance from the Li point to the 
center of the primary star has been divided into 10 equal 
segments in the center frame of Figure 1 , and the similar 
small circles on the Doppler tomograms are marked at 
velocity intervals corresponding to those distance inter- 
vals. Comparing the center and right frames of Figure 1, 
we see that the gas stream flow is truncated after travel- 
ing 6.5 small circles on the gas stream trajectory, which 
is equivalent to 65% of the distance from the Li point 
to the center of the primary star. The stream-star in- 
terac;tion occurs along the path of the gas stream right 
where the gas stream ends (at 65% of the distance from 
the Li point). However, Figure 1 (left) shows that the 
flow within the disk does not touch the star, since the 
gas within the disk (Feature 5) flows in the middle of the 
predicted locus of the disk. In addition, since the gas flow 
overlaps with the dashed circle in the tomogram, which 
corresponds to the Keplerian velocity at the Roche sur- 
face of the primary star, then the disk must extend out 
to the physical Roche surface of the star. The radius 
of the star, Rp = 3.0 Kq, q =0.29, the binary separa- 
tion, a = 1.23 X lO'^ km, and the Roche lobe radius, 
Rli = a(0.5 + 0.227 x loglO{q-^) = 11.0 Rq = 3.7 Rp. 
So, the inner edge of the accretion disk has a radius of 
~1.3i?p, while the outer edge of the disk has a radius of 
~3.7i?p. Both disk radii are measured from the center of 
the primary. 

A model of the U CrB binary is shown in Figure 9. 
This figure explains the various features detected in the 
3D tomograms of this binary. In the stream-like state 
(Fig. 9a), the geometry displays the case of direct im- 
pact without a disk and shows a luminous gas stream 
from the donor to the mass gainer. The 3D tomography 
results suggest that the stream impact site is displaced 
to the lower part of the star. In that case (see Fig. 9b) , 
it is reasonable that the splash and reflection of the gas 
stream (Feature 10 or FIO) can achieve moderate veloc- 
ity components {Vx ~ —60 km s^^, Vz ^ —200 km s~^; 
red lines). Beyond this splash site, the gas flows back to- 
wards the star (blue lines) and as a result. Feature 1 (Fl 
- annulus) has a positive displacement in ~ 80 km s~^ 
(see also Fig. 7 ). In the disk- like state (Fig. 9c), the 
geometry displays the case of direct impact with a disk 
but no intense luminous gas stream from the donor to 
the mass gainer (see also Figs. 2 and 6). The gas stream 
then interacts with the external part of the disk and the 
halo (gas located between the stars) and this interaction 
with the weak stream creates Feature 3 (F3) with veloc- 
ity components {V^ ^ —360 km s^^, ~ 400 km s""'^). 
However, the disk (and halo) decelerates or brakes and 
ejects gas with velocity components {Vx ~ 80 km s~^, 
Vz ^ 400 km s^^. The ejected gas can leave the system 
or fall back toward the star (and so it has negative Vz 
velocities). This is similar to the dynamical simulation 
results of Kuznetsov et al. (2001) (their Fig. 3a and 
Fig. 4). The remaining lower part of the gas stream 
strikes the disk but is too dilute to have a significant im- 
pact on the disk. This interaction is shown in Feature 



4 (F4), with velocity components {Vx ~ —580 km s~^, 
Vz ^ —200 km s~^). The Vz component for F4 is identi- 
cal to that for FIO, but the Vx component is 580 km 

because this represents the location where the stream 
ends, or perhaps because the stream impact occurs with 
the inner part of the disk, which has this value of Vx- 
The result is the formation of an annulus without any 
displacement in Vz- 

4. CONCLUSIONS 

The dynamical effects of gas flows perpendicular to 
the orbital plane and the vertical structure of the gas 
stream were considered in the classic theoretical work 
by Lubow & Shu (1976). Our 3D Doppler tomogra- 
phy results reveal that the process of mass transfer is a 
more complex phenomenon than expected. The accumu- 
lation of additional results from 3D Doppler tomography 
of different binaries will be useful in producing a better 
understanding of the multiple emission sources in these 
binaries and how the transverse flows influence the accre- 
tion process. The visualization of the features with the 
aid of 3D Doppler tomography represents a new step in 
data processing. It is especially useful in the identifica- 
tion and study of the transverse flows as well as features 
in the orbital plane. The comparison of the disk- like and 
stream-like states is illustrated in Figure 8 by means of 
a volume representation of the main central slices (V, = 
0) for 3D-1993 and 3D-1994 and also by means of color 
and contour images in this plane of the binary. It is clear 
that the disk-like state of U CrB in 1993 is distinct from 
the stream-like state in 1994. 

The main results of this work are (1) the application 
of a new tool to produce 3D Doppler tomograms, (2) 
the visualization of the disk-like and stream-like states 
of the Algol-type U CrB interacting binary, and (3) the 
discovery of gas flows moving across the orbital plane 
and the identification of the strongest emission sources 
(Features 1-10) in the 3D Doppler tomogram. 

1. The application of a new tool to produce 3D Doppler 
tomograms. The 3D Doppler tomograms created by us- 
ing the Radiastronomical Approach contain more infor- 
mation than the images produced with the standard 2D 
Filtered Back Projection Method. This new RA method 
represents a qualitative advance in data processing. The 
cascade of 2D images along the Vz direction have now 
provided a new understanding of the gas flows in the 
direction transverse to the orbital plane. This is true 
even without correcting for the artiflcial elongation of 
the sources in the Vz direction because of the lower res- 
olution in that direction. 

2. The visualization of the disk-like and stream-like states 
of the binary. The 2D Doppler tomograms of U CrB 
constructed from Ha line profiles collected in 1993 and 
1994 show two different states of the binary: a transient 
accretion disk representing the disk-like state (in epoch 
1993) and a bright gas stream flowing along predicted 
ballistic trajectory representing the stream-like state (in 
epoch 1994). These emission sources are even more dis- 
tinct in the 3D Doppler tomograms, and both states show 
evidence of a circumprimary emission source. Emission 
from the chromosphere of the cool secondary star is also 
visible especially in 1993 when the binary was in the disk- 
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like state. The multiple slices of the 3D tomogram show 
that the emission from the transient accretion disk is vis- 
ible in the Vz direction from Vz= -420 to +420 km s~^, 
although this range will be slightly smaller after the ad- 
justment for the lower resolution in the direction. The 
detection of negative 14 velocities as well as positive 
velocities associated with the disk emission suggests that 
the disk may be inclined to the direction of the orbital 
plane (the x,y) plane or we may be observing two arms 
of the disk, as observed in the cataclysmic variable IP 
Peg by Steeghs et al. (1997). 

The disk emission is not as strong as some of the other 
features. Moreover, the disk emission in the disk-like 
state is significantly weaker than the gas stream emis- 
sion in the stream- like state. These observations suggest 
that the disk is not as dense as the gas stream or the 
other compact features, and may explain why the tomo- 
gram of U CrB in its disk-like state does not show a 
symmetric disk which fills up its predicted locus in the 
{Vx,Vy) plane. The fact that the gas stream emission 
can dominate at certain epochs is further evidence that 
the gas stream is denser than the disk. These conclusions 
can be confirmed by comparing the hydrodynamic calcu- 
lations of Richards & Ratlifi^ (1998) for two Algol-type 
binaries, (3 Per and TT Hya, with those of Sawada et al. 
(1986); Spruit (1987); Kuznetsov et al. (2001) for the 
cataclysmic variable, IP Peg. Both TT Hya and IP Peg 
display dominant accretion disks over long timescales, 
while (3 Per shows no accretion disk. It is apparent that 
U CrB is different from IP Peg because the gas stream 
in U CrB is denser than the disk. This is why we can 
see the shock regions within the disk. However, in IP 
Peg, it seems that there is no gas stream, so the disk 
dominates the image so you can see the completed cir- 
cular pattern for the disk within the predicted locus of 
the disk. That is definitely not true for U CrB because 
the top part of the disk is missing, probably because the 
denser gas stream prevents the disk from completing the 
trip around the star. 

3. The discovery of transverse gas flmus and the identi- 
fication of the several emission sources. The reconstruc- 
tion of the 3D image of U CrB has clearly demonstrated 
that there exist strong gas flows which have high ve- 
locities. Therefore, we can conclude that these gas flows 
move in directions that are distinct from that represented 
by the standard 2D Doppler tomogram, which implicitly 
assumes that Vz = Q km s~^. In both states of this 
binary, a prominent source of emission is the circum- 
primary equatorial bulge on the primary star, which is 
produced when the gas stream strikes the surface of the 
star. In the disk-like state, the 3D image displays several 
additional compact emission sources. One of these is lo- 
cated along the gas stream trajectory, and is associated 
with the stream-disk interaction (14 ~ -300 km s~^, Vy ~ 
100 km s~i, Vz ^ 300 to 400 km s^^) and the other cor- 
responds to impact/splash region where the gas stream 
strikes the surface of the primary star {Vx ~ -580 km s~^, 
Vz ~ -300 to -hlOO km s'^). Such high velocity flows 
(jets) have been anticipated since the gas stream makes 
a high speed impact with the photosphere at speeds of up 
~500-600 km s~^, compared to the slower rotation of the 
mass gainer (Richards 1992; Blondin, Richards & Ma- 
linkowski 1995; Richards & Ratliff 1998). A second pair 



of bright features is found near the center of the (14, V^) 

plane at Vz velocities of ^ +200 to +500 km s""'^. These 
are found in a location similar to the Localized Region 
in /3 Per, RW Tau, and other direct-impact systems. 

In the stream-like state, there is evidence of a strong 
transverse flow with Vz velocity components of up to -300 
km s^^ (see Fig. 7). This gas jet may arise as a result 
of the stream-star interaction and overlaps with the lo- 
cation of this region in the 3D tomogram of U CrB when 
it was in the disk-like state. We have also found that 
the transverse flows associated with the circumprimary 
emission in the disk-like state flow in the opposite di- 
rection than when the binary is in the stream-like state. 
In the disk-Uke state, the flows have luminous Vz veloc- 
ity components in both positive and negative directions. 
Figure 6 shows that Features 3, 6, 7, and 8 have positive 
Vz c;omponents while Feature 4 has negative Vz veloci- 
ties. However, in the stream-like state, the transverse 
flow has a high negative Vz velocity of -300 km s~^ (see 
Feature 10 in Figure 7 for the 14= -60 km s~^ slice), 
while the portion with positive Vz velocities reaches a 
value of only ^100 km s"""^ (see Feature 1 in Figure 7 for 
the t4 = km s~^ slice). This suggests that the dense gas 
stream may penetrate through the rarified disk material 
and strike the star to form the circumbinary structure as 
an extension of the photosphere, but some of this mate- 
rial spreads out of the orbital plane, preferentially in the 
direction away from the observer, and at nearly seven 
times the synchronous velocity of the star. However, in 
the disk-like state, the disk material is denser than be- 
fore and so the gas stream does not have as strong an 
impact on the stellar surface. As a result, the gas is 
spread more evenly in both directions out of the orbital 
plane. Another explanation is that when the disk ma- 
terial circles the star and comes into contact with the 
denser gas stream, the interaction produces Features 6 
and 7 and displaces the stream-disk site at 3, on the 
outer part of the disk, so that the jet fiows in the op- 
posite direction. This explanation is consistent with the 
existence of an inclined accretion disk, and so we cannot 
rule out the possibility that the rotation axis of the pri- 
mary star may not be perpendicular to the orbital plane, 
as expected, and could be precessing. 

The calculation and analysis of 3D Doppler tomograms 
of additional systems will enable us to make a significant 
qualitative jump in the study of accretion in interact- 
ing binaries. The reconstruction procedure based on the 
Radioastronomical Approach will be influential in this 
future work. 
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Fig. 1. — 2D Doppler tomograms of U CrB-1993 in the disk-like state (left) and U CrB-1994 in the stream-like 
state (right) taken from Richards (2001), and a Cartesian model of the binary (center). 
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Fig. 2. — Slices of the 3D Doppler tomogram of U CrB from data collected in 1993. The images are shown for 
intensity levels: + 0.01, 0.05, 0.1, 0.3, 0.5, 0.7, 0.9 (red) and -0.01, -0.05, -0.1, -0.3 (blue), with velocities beyond the 
orbital plane of Vj = -420 km s^^ to = -f420 km s^^ in steps of 60 km s^^. 
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Fig. 3. — Four slices in the {Vx,Vy) plane for Vz = -240, -20 (maximum intensity), 240, 420 km s^^ showing the 
locations of the interesting features labeled 1-8. The levels of emission intensity in contour and color are the same as 
earlier in Figure 2. 
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Fig. 4. — The 3D Doppler tomogram of U CrB-1993 showing a cascade visuaHzation of the (V^jVy) shces from Vz 
= 240 to 540 km s~^ which display the interesting features labeled 3, 6, 7, 8. 
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Fig. 5. — The 3D Doppler tomogram of U CrB- 1994 showing a cascade visualization of the most interesting 
slices (adapted from Agafonov, Richards & Sharova (2006)). The resolution in the {Vx,Vy,Vz) directions correspond 
to velocities (30 x 30 x 110) km s~^. The three strongest Ha emission features are (1) the accretion annulus, (2) 
chromospheric emission, (9) the gas stream flowing along the ballistic trajectory from the LI point, and (10) the high 
velocity stream flowing in the Vz direction. 
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Fig. 6.— Cross sections in the {Vy,V^) plane in the 3D tomogram of U CrB -1993 for = -580, -360, -300, - 20, 
0, 80 km s"""^. The red contours correspond to the emission (positive intensities 0.05, 0.1, 0.3, 0.5, 0.7, 0.9) while the 
blue contours correspond to the absorption (negative intensities from -0.05 to -0.1). 
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Fig. 7. — Cross sections in the {Vy,Vz) plane in the 3D tomogram of U CrB-1994 for Vx = -60 and km s~^ 
adapted from Agafonov, Richards & Sharova (2006). The right frame illustrates emission from the equatorial bulge 
on the primary star. The left frame shows the transverse flow within the jet (10) as material is ejected from orbital 
plane away from the location where the gas stream strikes the surface of the mass gainer. The red contours correspond 
to the emission (positive intensities 0.05, 0.1, 0.3, 0.5, 0.7, 0.9) while the blue contours correspond to the absorption 
(negative intensities from -0.05 to -0.1). 
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Fig. 8. — Volume representation of the central slices (T4=0) and also their color+contour representation on the 
plane of the 3D Doppler tomogram calculated from the U CrB-1993 data compared to that found from the U CrB-1994 
data (from Agafonov, Richards & Sharova 2006). See intensity levels on the color scales. 
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Fig. 9. — A proposed model of the U CrB binary in Cartesian {x,z) coordinate space: a) the stream-hke state 
with the central direct stream impact on the mass gainer and without the disk; b) the displaced stream impact on the 
lower part of the mass gainer in the form of a jet (FIO; red line) and also gas falling onto the star (blue line); c) the 
disk-like state with a weak stream, which interacts with the disk. Features 3, 6, 7, 8 are the result of the diskstream 
interaction after the disk gas circles the star, but Feature 3 has a higher Vx velocity than the other features. The lower 
part of the stream interacts only with the inner part of the disk or with the photosphere of the star to produce feature 



F4, which is a jet with two velocity components {Vx ^ —580 km s ^ 
d) the same as in (c) but with the disk inclined. 



and Vz ^ 200 km s below the disk). Finally, 
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TABLE 1 

Characteristics and locations of prominent emission features in the 3D tomogram 



Location: central velocity or velocity range (km s ^) 



No. 


Emission Feature 


14 


Vy 




1 


Circumprimary emission 


(-60 to + 60) 


-60 (0 to -120) 


-20 (-180 to 180) 


2 


Chromospheric emission 


-80 to 80 


120 to 230 


-360 to 240 


3 


Stream-disk impact region 


-360 (-300 to -400) 


100 (0 to 200) 


420 (60 to 540) 


4 


Stream-star impact region 


-580 (-520 to -600) 


(-100 to 100) 


(-300 to +60) 


5 


Locus of the Accretion disk 


a 


a 


-420 to -f420 


6 


Localized Regions - Part 1 


-20 (-60 to +30) 


75 (30 to 120) 


200, 400 (150 to 540) 


7 


Localized Regions - Part 2 


75 (40 to 110) 


(-50 to -1-50) 


300 (200 to 540) 


8 


Feature near donor star 


80 (70 to 160) 


200 (130 to 270) 

b 


360 (200 to 540) 


9 


Gas stream flow 


b 




10 


High velocity transverse jet 


-60 





-400 to 



^corresponds to the predicted location of the feature on the 2D {Vx,Yy) plane, 
'^corresponds to the predicted ballistic trajectory of the gas stream. 



